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Exact Results for the Jarzynski Equality in Ising Spin Glass 
Models Derived by Using a Gauge Symmetry 

Chiaki Yamaguchi 
Kosugichou 1-359, Kawasaki 211-0063, Japan 



Exact results for the Jarzynski equality are derived for Ising spin glass models. The 
Jarzynski equality is an equality that connects the work in nonequilibrium and the differ- 
ence between free energies. The work is performed in switching an external parameter of 
the system. As the Ising spin glass models, the ±J model and the Gaussian model are in- 
jH , vestigated. For the ± J model and the Gaussian model, we derive exact lower bounds of the 

• exponentiated work for investigating the ferromagnetic phases and the multicritical points, 

and derive rigorous relations between the exponentiated work which have different quenched 
random configurations. For the Gaussian model, we derive the exact exponentiated work 
O^ ■ for investigating the spin glass phase. Exact results for the infinite-range models are also 

t^ ' obtained. The present results are obtained by using a gauge symmetry, and are related to 

points on the Nishimori lines which are special lines in the phase diagrams. The present 
(-H ' , results do not depend on any lattice shape, and a part of the present results instead depends 

pH ■ on the number of nearest-neighbor pairs in the whole system. 



§1. Introduction 



The theoretical studies of spin glasses have been widely done.'^''^' The spin glass 
models have the randomness and the frustration. The combination of the randomness 
and the effect of frustration causes various interesting dynamics as well as the static 
properties. For dynamical features, the aging phenomenon, the dynamical transition 
Q . for the distance between two spin configurations, and the problem of slow relaxation 

O \ are reported for example .l^™™ 

The Jarzynski equality is an equality that connects the work in nonequilibrium 

Cn . and the difference between free energies.'^ '1^ The work is performed in switching 

an external parameter of the system. The Jarzynski equality is also derived in the 



ff^ ■ Markov process with discrete time in Ref. 6), and it is pointed out in Ref. E]) that the 

f^ . Metropolis Monte Carlo method^' based on the Markov process with discrete time 

.\ is a suited example for applying the Jarzynski equality. Also, the Metropolis Monte 

^1j ' Carlo method and the related Monte Carlo methods are known as powerful methods 

^vq . to investigate the spin glass models on finite-dimensional lattices.'^ Therefore, the 

application of the Jarzynski equality and the related theories to the spin glass models 

may make its dynamical features clearer. Several exact results for applying the 

Jarzynski equality to the it J Ising spin glass model are also shown in Ref. |8]) , and it 

is suggested in Ref. fS) that the application of the Jarzynski equality and the related 

^_' theories to the spin glass models may be useful for avoiding the problem of slow 

relaxation, that reaching the equilibrium states is hard, in spin glass models. 

As the Ising spin glass models, the ± J modeP'^^ and the Gaussian modeP'''^''^ 
are investigated. These models have different quenched randomness. In Ref. [1]), for 
the lb J model, a number of studies by computer simulations are seen. In Ref. ED, for 
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the Gaussian model, a number of analytical studies are seen. The present results for 
the infinite-range models are also shown. 

In this article, a gauge transformation is used. There are special lines in the 
phase diagrams for several spin glass models, where the lines are called the Nishimori 
lines. ^-^ Several physical quantities and several bounds for physical quantities are 
exactly calculated on the Nishimori lines by using the gauge transformation.'^ The 
present results are related to points on the Nishimori lines. In the applications of 
many other methods, a lattice shape is supposed in advance, and the results are 
calculated on the lattice. On the other hand, in the method applied in this article, 
any lattice shape is not supposed in advance. The present results do not depend on 
any lattice shape, and a part of the present results instead depends on the number 
of nearest-neighbor pairs in the whole system. The present results are exact. 

In this article, we exactly analytically investigate the exponentiated work in the 
Jarzynski equality. The present results are useful for understanding the exact expo- 
nentiated work in the phase diagrams for the Ising spin glass models. In addition, 
the present results are useful for estimating the results by approximation approaches, 
and are useful for fixing bugs in computer programs for investigating the Ising spin 
glass models. 

This article is organized as follows. The Jarzynski equality and the models are 
explained in ^ Exact results for the itJ model are given in ^ and exact results 
for the Gaussian model are given in ^ This article is summarized in ^ 

§2. The Jarzynski equality and the models 

We explain the Jarzynski equality .'^''^''^ We define the inverse temperature of 
the reservoir as /3. Here, /3 = l/ksT, T is the temperature, and ks is the Boltzmann 
constant. We define a parameter for the strength of the exchange interaction between 
spins as J(> 0). We use a representation: K = f3J . By using Ka for K = Ka and 
-fCfe for K = Kfj, we consider a nonequilibrium process from Ka to K^. We assume 
that the initial and final states are in equilibrium, and the states in the process from 
Ka to Kh are in nonequilibrium. The Jarzynski equality is given by^'13'Eli 



-I3W = p-P'^F 



(2-1) 



where W is the work performed in the process from Ka to K^,, and the overbar 
indicates an ensemble average over all possible paths through phase space. AF is 
the difference between free energies given by AF = F(Kf,) — F{Ka), where F{Ka) 
is the free energy for K = Ka- The left-hand side of Eq. (|2-ip is the nonequilibrium 
measurements, and the right-hand side of Eq. ()2-ip is the equilibrium information. 
By using Eq. ()2Tp . one can extract the equilibrium information from the ensemble 
of nonequilibrium. Eq. ()2Tp does not depend on both the path from Ka to Xj,, 
and the rate at which the parameters are switched along the path. Note that, here, 
W is the work for changing J and a fixed /?. The work in switching an external 
parameter of the system in contact with a heat reservoir is originally supposed.'^''^ 
In Ref . Hj) , some particle-particle interactions which are turned on or off during the 
course of a molecular dynamics simulation are mentioned as a suited example for the 
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applications of the Jarzynski equality, and the present study may be exactly included 
in the example. As for the dynamics of the nonequilibrium process, the Markov 
process with discrete time is considered in this study, for exampleP Eq. (j2-l|) is 
equivalently given by 

^-mry •^■^' 

where Z[Ka) is the partition function for K = Ka- 

The Hamiltonian for Ising spin glass models, %, is given bj^^'Eli-llOJ 

n = -JY,^^,jS,Sj, (2-3) 

where (z,j) denotes nearest-neighbor pairs, Si is a state of the spin at the site i, and 
Si = ±1. Jtjj- is the strength of the exchange interaction between the spins at the 
sites i and j. Here, J > 0, and the value of J is switched from Ka/ P to Kb/p. The 
value of Tij is given with a distribution P{Tij). The ±J model and the Gaussian 
model are given by different distributions of Tij. 

The distribution P^^'^'{Tij) of tj.j for the it J model is given bj^''^ 

P^^'\r,,j) =VK,,1 + (1 - P) ^n,„-i , (2-4) 

where 5 is the Kronecker delta, p is the probability that the interaction is ferromag- 
netic, and 1 — p is the probability that the interaction is antiferromagnetic. By using 
Eq. ()2-4p . the distribution P^^-^\Tij) is rewritten as^ 



^K. 



(±J), 



P^^'\n,j)= " " Ij, , T,j = ±l, (2-5) 



2 cosh K^^^-^^ ' 



where Kp is given by 



K(-) = iln^. (2.6) 



lln^ 
2 1- 

When the value of Kp is consistent with the value of K, the line for K = Kp 
in the phase diagram is called the Nishimori line. The distribution P^'^'{Tij) of Tij 
for the Gaussian model is given bj^'l^'l^ 



^'°'<^''''^7S^'''^ 



2J2 



(2-7) 



and is rewritten as^ 



where Kp is given by 



Kf^ = 'f, . (2-9) 
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(G) (G) 

When the value of Kp is consistent with the value of K, the line for K = Kp in 

the phase diagram is also called the Nishimori line. 

A gauge transforniatioii^''^ given by 

Jnj — ^ JTijCTiaj , Si — )• SiCTi (2-10) 

is used, where cjj = ±1. The gauge transformation has no effect on thermodynamic 
quantities. ^^' By using the gauge transformation, the Hamiltonian 7i part becomes 
% ^%, and the distribution P^^'^'{Tij) part becomes 



A±J) ■ 
f:j^ *'^ 2^ (2 cosh K^^^ V^ 2^ (2 cosh if^^-^Vs ' 



^^^^^e^^^'E,,,)r.,,<..., ^ ^(^(±J)) 



where A^ is the number of sites, and Nb is the number of nearest-neighbor pairs in the 
whole system. By using the gauge transformation, the distribution pvGaussian)^^^ _^.^ 
part becomes 

(id) 

1 Y- / Y-^4- NBiK^Jif 



N 



B_ 






§3. Exact results for the ±J Ising spin glass model 

We derive exact results for the itJ Ising spin glass model. We multiply both 
sides of Eq. ([22]) by 



(2coshA:a)^s ' 
and perform the summation over {tjj}. Then, we obtain 

-^^ ^ / 2coshAA ^^ ^ e^°^^-^^->-'- Z{K,) 

^' J^'' V2coshKj ^^ f2coshKfe)^BZ(i^J' ^-^^^ 

where Wxa denotes the quenched random configuration average for Kp = K^- By 
performing the gauge transformation to the right-hand side of Eq. ()3Tp . we obtairP' 

2 cosh ATfe \ ^^ 



This equation is exact, is derived in Ref. [8]), and is called the Jarzynski equality for 
spin glass.l^ When Ka is fixed, the exponentiated work [e~^^]Ka is monotonically 
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Fig. 1. Two nonequilibrium processes and a schematic phase diagram for the ±J Ising spin 
glass model. Kp is a parameter for the concentrations of the ferromagnetic interaction and the 
antiferromagnetic interaction. K is a. parameter for the inverse temperature and the strength of 
the exchange interaction. The paramagnetic phase ('Para'), the ferromagnetic phase ('Ferro') 
and the spin glass phase ('Spin Glass') are depicted. The Nishimori line is also depicted as the 
dashed line. The nonequilibrium process for the exponentiated work [e~^^]Ka is depicted as 
the arrow 'A'. The nonequilibrium process for the exponentiated work [e~'^^]K^ is depicted as 
the arrow 'B'. 



increasing for Kf,. The exponentiated work [e~^^]Ka does not depend on any lattice 
shape, and instead depends on the number of nearest-neighbor pairs in the whole 
system, Nb- This equation is useful for estimating the results by approximation 
approaches, and is useful for fixing bugs in computer programs for investigating this 
model. 

When J — )• J l^fN and Nb = N{N — l)/2 are set, the model becomes the 
infinite-range model. We obtain 



-fiW] 



Ka 



exp 



{N-l){Kl-Kl] 



(3-3) 



for the infinite-range model with N ^ oo. This equation can be called the Jarzynski 
equality for the infinite-range itJ Ising spin glass model. When Ka is fixed, the 

^^]Ka with A^ — )• oo is monotonically increasing for Kb. 

we obtain [e~^^]Kp=o = 



exponentiated work [e " ]Ka with A^ 

When Ka = and Kb = K are set in Eqs. ([321) and ([3^ . 
exp[A'^B lii(coshi^)] for the finite-dimensional model and [e 



-fiWj 



K„=0 



exp[(A^ 



1)K /A] for the infinite-range model with N ^ oo. These result are results for 
p = 1/2, thus the randomness and the effect of frustration are largest in this model. 
In addition, these results are results for the nonequilibrium process from K = to 
K = K. Therefore, these results are results under simple conditions. 
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Fig.[T]shows two nonequilibrium processes and a schematic phase diagram for the 
lb J model. Kp is a parameter for the concentrations of the ferromagnetic interaction 
and the antiferromagnetic interaction, and is given in Eq. (|2-6|) . X is a parameter for 
the inverse temperature and the strength of the exchange interaction, and is given 
hy K = (3 J. The paramagnetic phase ('Para'), the ferromagnetic phase ('Ferro') 
and the spin glass phase ('Spin Glass') are depicted. The Nishimori line is also 
depicted as the dashed line. The nonequilibrium process for the exponentiated work 
[e~^^]Ka is depicted as the arrow 'A', and is mentioned above and in Ref. [8]). The 
nonequilibrium process for the exponentiated work [e~^'^]x(, is depicted as the arrow 
'B', and is mentioned below. 

The Cauchy-Schwarz inequality is given by 

[ E <{^i,^})yi{^^M' < E ^'({^..}) • E y'ii^i,^}) (3-4) 

in,]} {Ti,j} {n.j} 

for example, where x({rjj}) and y({Tjj}) are variables which depend on {tj^j}. When 
x{{Ti,j}) and y({rij}) are 






and 



y(K.}) = y(2coshK,)^BZ(i^„)' ('•') 

by applying the Cauchy-Schwarz inequality and the gauge transformation, we obtain 



2 cosh Ki, \ 



Nb 



2^^) SI-''"!-- P'^' 



The left-hand side of Eq. p-7p is a lower bound of the exponentiated work [e ^'^]ki 



The lower bound of the exponentiated work [e~P^]Ki does not depend on any lat- 
tice shape, and instead depends on the number of nearest-neighbor pairs in the 
whole system, Nb- This relation is rigorous, is useful for estimating the results by 
approximation approaches, and is useful for fixing bugs in computer programs for 
investigating this model. 

For the infinite-range model with iV — t- oo, we obtain 



exp 



{N-l){Kl-Kl) 



< [e-^'^U, . (3-8) 



The nonequilibrium process for the exponentiated work [e~^^]Ka is limited to 
the under left side of the Nishimori line in the phase diagram shown as the arrow 'A' 
in Fig. [TJ Therefore, Eqs. (j3-2|) and (j3-3p are useful for investigating the spin glass 
phase. The nonequilibrium process for the exponentiated work [e~^^]Ki, is limited 
to the upper right side of the Nishimori line in the phase diagram shown as the arrow 
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'B' in Fig. [TJ Therefore, the left-hand sides of Ineqs. ()3-7p and (|3-8p are useful as 
lower bounds in investigating the ferromagnetic phase and the multicritical point for 
the paramagnetic phase, the ferromagnetic phase and the spin glass phase. 
By using Eqs. (|3-2p and ()3-7p . we obtain 



[e-nK. < [e-^'^W ■ 



(3-9) 



The relation (|3-9p does not depend on any lattice shape, and also holds for the 
infinite-range model. This relation is a rigorous relation between the exponentiated 
work which have different quenched random configurations. This relation is useful 
for understanding the exact exponentiated work in the phase diagram. This relation 
exactly shows that the exponentiated work of the nonequilibrium process shown as 
the arrow 'B' in Fig. [T]is greater than or equal to that of the nonequilibrium process 
shown as the arrow 'A' in Fig. [TJ 



§4. Exact results for the Gaussian Ising spin glass model 



We derive exact results for the Gaussian Ising spin glass model. We multiply 
both sides of Eq. ^T^ by 

^2 



(2vrJf: 



exp 



^2Jl 2 ^^^^ 

{id) ^ (id) 



«J 



and perform the integration over {tij}. Then, we obtain 



-pWj 



Ka 



(2vrJf)- 



e 2 



(II^^mO 



■E(ij> :^+Ko.T,{i,j)-ri,i 



i^d) 



Z{Ka 



(4-1) 



By performing the gauge transformation to the right-hand side of Eq. ()4-ip . we obtain 

■Nb{kI-kI)jI^ 



[e-P^]K^ = exp 



(4-2) 



This equation is exact, and can be called the Jarzynski equality for the Gaussian 
Ising spin glass model. When Ka is fixed, the exponentiated work \e~^'^]Ka is mono- 
tonically increasing for i^^. The exponentiated work [e~^'^]x^ does not depend on 
any lattice shape, and instead depends on the number of nearest-neighbor pairs in 
the whole system, Nb- This result is useful for estimating the results by approxima- 
tion approaches, and is useful for fixing bugs in computer programs for investigating 
this model. 

When J — )■ J/yN and Nb = N{N — l)/2 are set, the model becomes the 
infinite-range model, i.e., the Sherrington-Kirkpatrick model.'^ We obtain 

■{N-l){Kl-Kl)Jp 



[e-P^]K^ = exp 



(4-3) 
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Fig. 2. Two nonequilibrium processes and a schematic phase diagram for the Gaussian Ising spin 
glass model. Kp is a parameter for the concentrations of the ferromagnetic interaction and the 
antiferromagnetic interaction. K is a. parameter for the inverse temperature and the strength of 
the exchange interaction. The paramagnetic phase ('Para'), the ferromagnetic phase ('Ferro') 
and the spin glass phase ('Spin Glass') are depicted. The Nishimori line is also depicted as the 
dashed line. The nonequilibrium process for the exponentiated work [e~^^]_ff„ is depicted as 
the arrow 'A'. The nonequilibrium process for the exponentiated work [e"''^]^^ is depicted as 
the arrow 'B'. 



for the Sherrington-Kirkpatrick model. This equation can be called the Jarzynski 
equality for the Sherrington-Kirkpatrick model. When Ka is fixed, the exponentiated 
work [e~^'^]Ka is monotonically increasing for Kb. 

When Ka = and Kb = K are set in Eqs. (021) and g31), we obtain [e~>^^]Kp=o = 
exp[NB{KJi)'^/2] for the finite-dimensional model and [e~l^^]Kp=o = exp[(A^ — 
l)(i^Ji)^/4] for the infinite-range model. These result are results for Jq = 0, thus the 
randomness and the effect of frustration are largest in this model. In addition, these 
results are results for the nonequilibrium process from K = to K = K. Therefore, 
these results are results under simple conditions. 

Fig. [2] shows two nonequilibrium processes and a schematic phase diagram for 
the Gaussian model. Kp is a parameter for the concentrations of the ferromagnetic 
interaction and the antiferromagnetic interaction, and is given in Eq. ()2-9p . K is a. 
parameter for the inverse temperature and the strength of the exchange interaction, 
and is given hy K = (3J. The paramagnetic phase ('Para'), the ferromagnetic phase 
('Ferro') and the spin glass phase ('Spin Glass') are depicted. The Nishimori line is 
also depicted as the dashed line. The nonequilibrium process for the exponentiated 
work [e~l^'^]Ka is depicted as the arrow 'A', and is mentioned above. The nonequi- 
librium process for the exponentiated work [e~^^]Kt is depicted as the arrow 'B', 
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and is mentioned below. 

The Cauchy-Schwarz inequality is given by 



oo 
oo 



Yl dTij) x{{Tij})y{{Tij})f 



(»J) 



/OO /"OO /"OO /"OO 

•••/ (n^^M)^'(K^»- / •••/ (11 ^^^^^^^'(Ki}) (4-4) 
-OO J— oo /. .\ J —oo J —oo I- -v 



for example. When variables 2;({tjj}) and y({Tjj}) are 



a^({T-i,j}) 



1 






N (2vrJ?)' 



(4-5) 



and 



y({-^ij}) 



~ S<i.i> 57f ^— T"^— +-f^« E(i,j) ^i,j 2'(Kb) 



(4-6) 



by applying the Cauchy-Schwarz inequality and the gauge transformation, we obtain 



exp 



^b{KI 



KlVl 



< [e-P^] 



Kh 



(4-7) 



The left-hand side of Eq. (|4-7p is a lower bound of the exponentiated work [e ^'^]ki,- 
The lower bound of the exponentiated work [e"^^];^-^ does not depend on any lat- 
tice shape, and instead depends on the number of nearest-neighbor pairs in the 
whole system, Nb- This relation is rigorous, is useful for estimating the results by 
approximation approaches, and is useful for fixing bugs in computer programs for 
investigating this model. 

For the Sherrington-Kirkpatrick model, we obtain 



exp 



iN-l){Ki 



Kl)Jl 



< [e-^^]x. • 



(4-8) 



The nonequilibrium process for the exponentiated work [e~l^^]Ka is limited to 
the under left side of the Nishimori line in the phase diagram shown as the arrow 'A' 
in Fig. [2j Therefore, Eqs. (j4-2p and (j4-3p are useful for investigating the spin glass 
phase. The nonequilibrium process for the exponentiated work [e~^'^]i<-(, is limited 
to the upper right side of the Nishimori line in the phase diagram shown as the arrow 
'B' in Fig. [2l Therefore, the left-hand sides of Ineqs. ()4-7p and ()4-8p are useful as 
lower bounds in investigating the ferromagnetic phase and the multicritical point for 
the paramagnetic phase, the ferromagnetic phase and the spin glass phase. 

By using Eqs. ()4-2p and ()4-7p . we obtain 



[e-P'^U. < [e-/^^]^, 



(4-9) 



10 C. Yamaguchi 

The relation ()4-9p does not depend on any lattice shape, and also holds for the 
Sherrington-Kirkpatrick model. This relation is a rigorous relation between the expo- 
nentiated work which have different quenched random configurations. This relation 
is useful for understanding the exact exponentiated work in the phase diagram. This 
relation exactly shows that the exponentiated work of the nonequilibrium process 
shown as the arrow 'B' in Fig. [2]is greater than or equal to that of the nonequilibrium 
process shown as the arrow 'A' in Fig. [2j 

§5. Summary 

Exact results for the Jarzynski equality were derived for Ising spin glass models. 
As the Ising spin glass models, the it J model and the Gaussian model were investi- 
gated. For the itJ model and the Gaussian model, we derived exact lower bounds 
of the exponentiated work for investigating the ferromagnetic phases and the multi- 
critical points, and derived rigorous relations between the exponentiated work which 
have different quenched random configurations. For the Gaussian model, we derived 
the exact exponentiated work for investigating the spin glass phase. Exact results 
for the infinite-range models were also obtained. The present results were obtained 
by using a gauge symmetry. The present results are related to points on the Nishi- 
mori lines. The present results do not depend on any lattice shape, and a part of 
the present results instead depends on the number of nearest-neighbor pairs in the 
whole system. 

The present results are useful for understanding the exact exponentiated work 
in the phase diagrams for the Ising spin glass models. In addition, the present results 
are useful for estimating the results by approximation approaches, and are useful for 
fixing bugs in computer programs for investigating the Ising spin glass models. 

References 

1) N. Kawashima and H. Rieger, "Recent Progress in Spin Glasses", in Frustrated Spm Sys- 
tems, ed. H. T. Diep (2004). 

2) M. Mezard, G. Parisi and M. A. Virasoro, Spin Glass Theory and Beyond (World Scientic, 
Singapore 1987). 

3) B. Derrida and G. Weisbuch, Eurphys. Lett. 4 (1987), 657. 

4) C. Jarzynski, Phys. Rev. Lett. 78 (1997), 2690. 

5) C. Jarzynski, Phys. Rev. E 56 (1997), 5018. 

6) G. E. Crooks, J. Stat. Phys. 90 (1998), 1481. 

7) N. Metropolis, A. Rosenbluth, M. Rosenbluth, A. Teller and E. Teller, J. Chem. Phys. 21 
(1953), 1087. 

8) M. Ohzeki and H. Nishimori, J. Phys. Soc. Jpn. 79 (2010), 084003. 

9) H. Nishimori, J. Phys. C 13 (1980), 4071; Prog. Theor. Phys. 66 (1981), 1169. 

10) D. Sherrington and S. Kirkpatrick, Phys. Rev. Lett. 35 (1975), 1792. 

11) G. Toulouse, Commun. Phys. 2 (1977), 115. 



